Renormalization persistency of tensor force in nuclei 



Naofumi Tsunoda 

Department of physics, the University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, Japan 

Takaharu Otsuka 
Department of physics and Center for Nuclear Study, 
the University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, Japan 
National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, MI, 48824, USA 

Koshiroh Tsukiyama 
Department of physics and Center for Nuclear Study, 
the University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, Japan 



o 

(N 
D 
00 



(N 
> 



Morten Hjorth- Jensen 

Department of Physics and Center of Mathematics for Applications, University of Oslo, N-0316 Oslo, 
National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, MI, 



Norway 
, USA 



In this work we analyze the tensor-force component of effective interactions appropriate for nuclear 
shell-model studies, with particular emphasis on the monopole term of the interactions. Standard 
nucleon-nucleon (NN) interactions such as AV8' and xN^LO are tailored to shell-model studies 
by employing Viowk techniques to handle the short-range repulsion of the NN interactions and by 
applying many-body perturbation theory to incorporate in-medium effects. We show, via numerical 
studies of effective interactions for the sd and the pf shells, that the tensor-force contribution to the 
monopole term of the effective interaction is barely changed by these renormalization procedures, 
resulting in almost the same monopole term as the one of the bare A*' A'' interactions. We propose to 
call this feature Renormalization Persistency of the tensor force, as it is a remarkable property of 
the renormalization and should have many interesting consequences in nuclear systems. For higher 
multipole terms, this feature is maintained to a somewhat smaller extent. We present general 
intuitive explanations for the Renormalization Persistency of the tensor force, as well as analyses 
of core-polarization terms in perturbation theory. The central force does not exhibit a similar 
Renormalization Persistency. 

PACS numbers: 21.30.Fe, 21.60.Cs 
Keywords: tensor force, effective interaction 
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I. INTRODUCTION 

The nucleon-nucleon (NN) interaction is normally 
modelled in terms of several components, such as a cen- 
tral force, a spin-orbit force or a tensor force. These 
mathematical terms accommodate our phenomenological 
knowledge of the strong interaction, vifhich, when used 
in a nuclear many-body context is subjected to different 
renormalization procedures. For the nuclear many-body 
problem, a given renormalization procedure leads to the 
derivation of an effective interaction, starting from a bare 
realistic NN interaction. The so-called bare NN interac- 
tions exhibit a strong coupling between low-momentum 
and high-momentum degrees of freedom generated from 
short-range details of the interaction. By "bare" we 
mean that the abovementioned strong coupling is left 
untouched. This coupling is included only implicitly, via 
various renormalization procedures, in the effective inter- 
actions used in for example shell-model studies. 

As an example of bare NN interactions, the Argonne 
interactions (AV), which are defined in terms of local op- 
erators in coordinate space, show a strong short-range 
repulsion [l|, [2] • The resulting strong coupling between 



low- and high-momentum modes makes the many-body 
problem highly non-perturbative. On the other hand, 
in shell-model calculations, the employed effective in- 
teractions are defined for a specific configuration space 
(a strongly reduced Hilbert space), normally called the 
model space. Therefore, the effective interactions for the 
shell model should be renormalizcd to include the effects 
of virtual excitations to the configurations not included 
in the model space. 

Although the properties and the effects of the full inter- 
action and various renormalized interactions have been 
investigated extensively over the years, we feel that there 
are still important features of the nuclear interaction 
which deserve some special attention. In particular, we 
show here via several numerical studies, that the tensor 
force component of the bare nuclear interaction is left al- 
most unaffected by various renormalization procedures. 
The monopole component of the tensor force, a compo- 
nent of great interest in studies of shell evolution (see 
discussion below) in nuclei toward the drip lines, is left al- 
most unchanged under various renormalizations. This al- 
lows us thereby to extract simple physics interpretations 
from complicated many-body systems. In this work we 
label such a lack of renormalization influence as Renor- 
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malization Persistency, in short just R-Persistency. The 
R-Persistency is a property exhibited by specific terms 
of the original nuclear Hamiltonian that arc not affected, 
or barely affected, by the renormalization procedure. 

On the experimental side, present and future radioac- 
tive ion-bcam facilities have made it possible to perform 
experiments that explore nuclei far from the stability line 
of the nuclear chart. Many unexpected and new phenom- 
ena have been observed in such experiments carried out 
at radioactive ion-beam facilities worldwide. One of the 
most striking results is the breaking of the conventional 
shell structures in neutron-rich nuclei. Such shell evolu- 
tion, unexpected in the past, is known by now to occur 
mainly due to an unbalanced neutron to proton ratio 
and specific rbital-dependent components of the nuclear 
forces. In particular, the nuclear tensor force plays a key 
role here, as proposed by one of the authors [3, 3- One 
of the most useful quantities to probe the effect of ten- 
sor force is the so-called monopole matrix element. The 
monopole matrix element Q of the two-body interac- 
tion between two single-particle states labelled j and j' 
and total two-particle isospin T is defined as 

yT _ j:A2J+l){jf\V\jf)jT . . 

E,7(2J + 1) ' ^' 

for j ^ j' [13]. Here {■■\V\--)jt denotes the anti- 
symmetrized two-body matrix element coupled to total 
angular momentum J and total isospin T. The monopole 
matrix element is crucial for shell evolution, because it 
affects the effective single particle energy linearly. For in- 
stance, if neutrons occupy the single-particle state 
j', they shift the effective single particle energy of protons 
in the state j as follows, 

^^pW = ^(^y=° + ^y=') (2) 

where Aep(j) represents the change of the effective sin- 
gle particle energy of protons in the single-particle state 
j. When we consider the tensor- force contribution, the 
monopole matrix elements always have different signs be- 
tween a pair of spin-orbit partners. For example, the in- 
teraction matrix elements Vj-^j' and Vj_^ji have opposite 
sign. Here we define j> and j< to represent the spin-orbit 
partners, that is, j> =1 + 1/2 and j< = I ~ 1/2, where 
I stands for the orbital momentum of a given single- 
particle state. In this case, the tensor-force changes the 
spin-orbit splitting between j> and . The shell struc- 
ture is also altered, in particular if we have a sizable 
number of neutrons in the single-particle state j' . 

In previous studies Q, the tensor-force component in 
effective interactions for shell-model calculations was, for 
the sake of simplicity, modelled via the exchange of tt and 
p mesons only. To a large extent, this yields results close 
to the tensor force in realistic NN interactions. In fact, 
this model describes rather well the experimental data in 
several mass regions However, it is far from trivial 
that the tensor force in effective interactions for the shell- 



model can be considered to be given by the exchange of 
TT and p mesons only. 

The aim of this article is thus to investigate the R- 
Persistency of the nuclear tensor force and understand 
the validity of the above assumption through theoretical 
studies, based on realistic NN interactions and micro- 
scopic theories for deriving effective interactions, focus- 
ing on the effective interaction for the shell model j2]| . 

This work is organized as follows. First we briefly re- 
view the theory for constructing effective interactions in 
Sec. HH In Sec. Iml and Sec. El the R-Persistency of 
the monopole part of the tensor force from various ap- 
proaches to the effective interactions will be discussed. In 
Sec. IVII we present not only the monopole part but also 
the two-body matrix elements including multipolc part 
of the various effective interactions and discuss their ten- 
sor force components. For the sake of completeness, we 
include analyses using other NN interactions in Sec. I VIII 
The last section contains our conclusions. 



II. CONSTRUCTION OF THE EFFECTIVE 
INTERACTION FOR THE SHELL MODEL 

The aim of this section is to give a brief sketch of the 
theoretical methods we employ in our analyses of the nu- 
clear force. To construct the effective interactions for 
the nuclear shell model, we use many-body perturba- 
tion theory (MBPT). However, as inputs to MBPT, we 
cannot use bare realistic NN interactions directly, since 
their high-momentum components make MBPT non- 
convergent 01- We integrate out this high- momentum 
components employing a renormalized interaction de- 
fined only in the low-momentum space below a certain 
sharp cutoff A and designed not to change two-body ob- 
servables like NN scattering data Q. This recipe defines 
a cutoff dependent family of interactions, normally la- 
belled as T^owfc(A), which to be more specific, can be 
written as 

Mowfc(A) = PAl/barc^A + 6Vct{A), (3) 

where Pa indicates a projection operator onto the low- 
momentum space below A. The term SVct{A) represents 
the correction term coming from the renormalization pro- 
cedure. In other words, PA^baroPA is a simple projection 
to a low- momentum space, while 6Vct{A) emerges as a 
result of the chosen renormalization procedure. By con- 
struction, T^owfc(A) approaches the original NN interac- 
tion in the limit A — > oo. A complete renormalization 
scheme would generate higher-body forces as well, such 
as three-body and four-body forces, VJjn and V4N, respec- 
tively. In this work we limit ourselves to two-body (VJjn) 
interactions only. Thus the cutoff dependence of phys- 
ical quantities can be used to assess the error made by 
omitting more complicated many-body forces. The term 
PAVharcPA should coutain the long-range part of one-pion 
exchange interaction as a major component. 
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Next, we proceed to MBPT. The low-momentum in- 
teraction Viovjk is a good starting point for MBPT be- 
cause we can avoid the difficulty caused by the strong 
short-range repulsion. For a degenerate model space, the 
effective interaction VeS can be written as 



decomposition employed in for example Refs. |12h1 



(4) 



where Q{Eo) is the so-called Q-hox, defined as 



Q{Eo) = PHiP + PHiQ 



1 



Eo - QHQ 



QHiP. (5) 



Here, the Hamiltonian is divided into an unperturbed 
part Hq and an interaction part Hi , H = Hq + Hi and 
the model space is set to be degenerate with respect to 
the unperturbed Hamiltonian Hq with energy Eq. The 
integration symbols in Eq. (|4]) represent the inclusion to 
infinite order of so-called folded diagrams, see Refs. [olIToj 
for details. The Q-box is given by diagrams which are 
valence-linked and irreducible, while Q' indicates that 
only diagrams which are of second- or higher-order in 
terms of the interaction Hi , arc included, sec for example 
Fig.ffl 

We can solve Eq. (|4]) by the following iterative formula 



K^;) = Q{Eo) + V Q,n{Eo){V,T"} 



E 

m— 1 



(6) 



where Q^(Eo) ^ ^A^^"^^ ' 
ml du 



In this work, we take 



into account diagrams up to second or third order in the 
interaction Hi for the calculation of the Q-box. 




(c) 



FIG. 1: Examples of diagrams to second-order in the inter- 
action Hi included in the Q-box. The diagrams are referred 
to as (a) particle-particle ladder, (b) hole-hole ladder and (c) 
core-polarization, respectively. 



By using this two-step method, we can start from an 
arbitrary bare realistic NN interaction. We calculate ef- 
fective interactions starting from AV8' P, and the chi- 
ral xN'^LO interaction Results using the AV8' in- 

teraction are shown in the following sections while our 
results obtained with the xN'^LO interaction are shown 
in Sec. IVIII for the sake of completeness. 

Finally, to extract the tensor component from the ob- 
tained effective interactions, we employ the spin-tensor 



{abLS\Vp\cdL' S'),,,T = (-l)'^'p 



L S J' 
S' L' p 



(7) 



where ( • • LS*]!^! • • L'S')jt denotes the LS'-coupled ma- 
trix element of the effective interaction. Here a (as well as 
bed) is shorthand for the set of quantum numbers (ria, la)^ 
etc. The operator Vp is defined as the scalar product 
Vp = U^P^ ■ X^P\ where U'^p^ and X^p^ are irreducible 
tensors of rank p, applying to operators in both spin and 
coordinate space. The tensor component is extracted by 
setting p = 2 in Eq. ([7]). Finally, in the above equation 
we have defined p = 2p+ 1 and J = 2 J + 1. 



III. TENSOR FORCE IN LOW-MOMENTUM 
INTERACTION V^iowfe 

We now present results obtained by the theoretical 
methods described in the previous section. Figure [2] 
shows the monopole part of the tensor-force of the renor- 
malized Viowfe interaction derived from the Argonne V8' 
(AV8') potential for the sd-shell and the p/-shcll. The 



cutoff value A varies from 1.0 fm"' to 5.0 fm' 



Here, 



we employ units where c = h = h?/m — 1. The typical 
value of the cutoff is determined by the best reproduc- 
tion of the binding energies of ^H and ''He. The result- 
ing cutoff value lies around 2.0 fm~^ [l^. A too small 
cutoff A (for example 1.0 fm~^ in momentum space) can- 
not resolve the necessary degrees of freedom. Since the 
Compton length of the pion is approximately 0.7 fm, a 
cutoff A = 1.0 fm~^, which corresponds to 1.0 fm in co- 
ordinate space, is too small to resolve the exchange of 
a pion. Although the resulting renormalized interaction 
Viovjk with A = 1.0 fm~^ may not contain an appropri- 
ate tensor force for shell- model calculations, wc include 
its result in Fig. [2] and subsequent similar figures for the 
sake of completeness. 

We now present the results for the cutoff values A = 
1.0,2.1 and 5.0 fm^^ in Fig. [21 The matrix elements are 
calculated using a harmonic oscillator basis with Huj = lA 
and 11 MeV for the sd-shell and the p/-shell, respec- 
tively. Except for a very low (and thereby unreasonable) 
cutoff value A = 1.0 fm~^, one finds, both in the sd-shell 
and the p/-shell, that the monopole part of the tensor 
force of Viowfc has almost no cutoff dependence, and has 
almost the same strength as that of the original NN in- 
teraction. Thus, within the usual values of the cutoff, 
we can sec that the monopole part of the tensor force 
fulfills the R-Persistency almost perfectly with respect to 
the renormalization of the short-range part of the NN 
interaction. 

Wc look now into the robustness and the generality 
of the features discussed above. For this purpose, we 
consider the relative motion of two interacting nucleons. 
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FIG. 2: (color online.) Tensor-force monopole component of low-momentum interaction Viowfc as function of the cutoff parameter 
A for (a.) T — forces in p/-shell, (b) T = 1 forces in p/-shell, (c) T — forces in sd-shell and (d) T — 1 forces in sd-shell. The 
cutoff parameter A of Viowfe varies from 1.0 fm"^ to 5.0 fm"\ 



The orbital angular momentum of the relative motion 
can be L=0 (S), 1 (P), 2 {D), etc. If the tensor force is 
acting between two states, there is no coupling between 
two S states, because the relative coordinate operator in 
the tensor force is of rank 2. The S'-to-S' coupling is thus 
zero. This results in strongly suppressed contributions to 
the tensor force from the short-range part of the relative- 
motion wave function, since a good fraction of the short- 
range repulsion stems from S waves. Partial waves higher 
than S waves carry also a centrifugal barrier component 
which results in smaller short-range contribution to the 
tensor force relative to S waves. Thus, changes of the po- 
tential at short distances do not affect matrix elements of 
the tensor force for low momentum states. This seems to 
be the basic reason why the tensor force remains almost 
the same throughout the renormalization procedure. In 
other words, there is a sound reason to expect the R- 
Persistcncy for the tensor force regarding the treatment 
of the short-range correlation. On the other hand, the 
present argument may not be applied to other parts of 
the nuclear force such as the central force. 

The second term SV of Eq. ([3]) is due to the renormal- 
ization. It includes for example the central-force com- 
ponent at intermediate inter-nucleon distances, and may 
affect, in principle, the tensor force as well. The first 
term, PAVbaic-PA, is equal to the bare tensor force in 
the limit of A — >■ oo by definition. In this limit SV is 



zero. Since matrix elements of the tensor force, partic- 
ularly for low-momentum states, are not affected much 
by the short-range modification, the effect of the tensor- 
force component in the first term of Eq. ([3]) remains the 
same to a large extent, even with finite A values, un- 
less it becomes extremely small. The fact that the R- 
Persistency is almost fulfilled in numerical calculations 
(as we can see in Fig. [2]) implies therefore that the sec- 
ond term 6V results in small contributions to the tensor 
force, or does not change the long-range part of the ten- 
sor force. The origin of the weak tensor force component 
in SV can be understood by the arguments presented in 
Sec. |Vl arguments which are based on the close relation 
between the Mowfc renormalization process and contribu- 
tions from MBPT that represent long-range corrections, 
as discussed in d, [l^ as well. We shall come back to this 
point in Sec. |Vl 



IV. RENORIMALIZATION OF THE CENTRAL 
FORCE 

Contrary to the tensor force, it can be seen from our 
numerical studies that the central force does not fulfill 
the R-Pcrsistcncy, and is indeed affected strongly by the 
renormalization procedure due to the short-range part of 
the NN interaction. This is reflected in a much stronger 
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cutoff dependence as well. The central-force monopole 
part of SV in Eq. ^ is thus not smaU. 
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FIG. 3: (color online.) Central-force component of the 
monopole term of the bare AV8' (bare in the figure), Viowfe 
(TS) and Viowfc (full) for the sd-shell, see text for further de- 
tails and discussions. The central-force component is obtained 
using the decomposition of Eq. 0. The effect of the renor- 
malization on the short-range tensor force is also shown. The 
cutoff value is chosen as A — 2.1 fm~^. 



Figure [3] shows the monopole part of the central force 
of the bare AV8' potential obtained by the decomposition 
of Eq. ([7]). In this figure we show also the corresponding 
central-force monopole component using the Viowfc renor- 
malized interaction originating from the AV8' potential, 
labelled by "full" in the above figure. We show also re- 
sults where the tensor force component has been sub- 
tracted from the bare NN interaction in the renormaliza- 
tion procedure, labelled by "TS" in the figure. What we 
can see in Fig.[3]is the effect of the renormalization due to 
the short-range part of the bare realistic NN interaction. 
The difference between "bare AV8' " and "Mowfe (TS)" 
lies mainly in the renormalization due to the short-range 
part of the central force, as the tensor force is subtracted 
in "Mowfc (TS)". On the other hand, the difference be- 
tween "Viowfe (TS)" and "Vowfc (full)" comes solely from 
the renormalization due to the short-range part of the 
tensor force. 

In the T = channel, the effect of the renormalization 
procedure on the short-range part of the tensor force is 
comparable to that of the central force, while in the T = 1 
channel this effect is almost negligible. This is a quite 
remarkable feature. Let us discuss this feature in some 
detail by considering the Schrodinger equation for the 
deuteron. The deuteron has isospin T = 0, spin S = 1, 
orbital momentum L ~ (S'-wave) and total angular 
momentum J = 1. There is a small admixture of D- 
waves as well, leading to the following coupled differential 



equations for the deuteron 



M dr2 
Yl dr2 



Vcu{r) 

6h^ 



+V8VTu{r) = Edw{r 



Vc - 2Vt - 3Vls ) w{r) 



(8) 



where u(r) and w{r) are the radial wavefunctions of the 
5-wave and the D-wave, respectively. The potentials 
Vc , Vls and Vr are the central, spin-orbit and tensor 
forces, respectively. Knowing the solution of Eq. 
we can integrate out the Z?-wave degrees of freedom and 
obtain the following effective central force 



Veff (r; 35i) = Vc(r; ^^i) + AVeg (r; ^^i), 

AVeff(r;35i) = y8VT(r)44' 

u[r) 



(9) 



The effective central force AV-ff is comparable to Vc in 
strength and it makes the ^Si channel the most attrac- 
tive one [13, [13 ■ This effective central force makes the 
deuteron bound for the '^Si channel. We can regard this 
equation as a special case of Eq. (jS]). The effective cen- 
tral force comes from a second-order effect due to tensor 
force, since both the initial and the final state have or- 
bital angular momentum 0. As a consequence, the effec- 
tive interaction for the T = channel is enhanced by the 
renormalization procedure due to the short-range part of 
the tensor force. This is however not the case in T = 1 
channel. It reflects the property of the deuteron, which 
is the only bound two-nucleon system. A similar mech- 
anism may also explain the strong cutoff dependence of 
the Vowfc interaction seen in the T = channel. 



V. TENSOR FORCE IN EFFECTIVE 
INTERACTION FOR THE SHELL MODEL 

We discuss here the tensor-force component in the ef- 
fective interactions for the shell model, using the decom- 
position of Eq. ([T)). We have calculated effective inter- 
actions for the shell model (V^^^) using many-body per- 
turbation theory (MBPT) by considering the Q-box up 
to second and third order with folded diagrams included 
as well, starting from a renormalized Vowfc interaction. 
The cutoff value used in the Vowfc calculation is set to 
A = 2.1 fm~^ [l^. The model space (P-space) is cho- 
sen to be the full sd-shell or the full p/-shell. In the 
construction of the Vowfe interaction, we renormalize the 
strong short-range repulsion of the NN interaction, and 
in MBPT we include further effects of truncations of the 
model space. The Q-box is calculated by considering 
valence-linked and connected diagrams with unperturbed 
single particle energies of the harmonic oscillator. The 
oscillator energy hut is set to be 14 MeV and 11 MeV 
for the sd-shell and the p/-shell effective interactions, re- 
spectively. Degenerate perturbation theory is employed 
in constructing the effective interactions. 
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of 



Since the Q-space is defined as the complement of 
the P-space, intermediate states arising in each diagram 
should be taken up to infinitely high oscillator shells. In 
our case, using a low-momentum interaction VIowfe with 
A = 2.1 fm^^, full convergence of the monopole part of 
VcS is obtained with approximately 8—10 hio excitations 
in each diagram which makes up the Q-box. 

Figure 2] shows the monopole part of the tensor force 
defined for the sd-shell or the p/-shell. As a 
general trend, one can see again that the monopole part 
of the tensor force of V^f^ fulfills our R-Persistency hy- 
pothesis to a good extent both in the sd-shell and in the 
p/-shell. Since the first order Q-box is just the Mowfc in- 
teraction, the results mean that the monopole part of the 
tensor force is dominated by the first-order term in the 
Q-box and the contributions from second or higher-order 
terms are remarkably small. These results can be un- 
derstood by considering the specific angular momentum 
structure of the tensor force, which is a scalar product 
of two rank 2 tensors in spin and coordinate spaces. In 
a perturbative correction to second or higher-order, such 
a complicated structure is smeared out and the resulting 
interaction consists mainly of a central force contribu- 
tion. Therefore, as for the tensor-force component in the 
monopole interaction, it is the first-order contribution 
which is the dominant one. 

To elucidate why higher-order terms in many-body 
perturbation theory are small, we consider as an example 
a contribution from second order in the interaction, by 
far the largest higher-order term. 

The Hamiltonian causing the present second-order per- 
turbation can be written as 



ip)) 



(10) 



p=0,l,2 



where Wp represents an interaction strength, C/^^-' and 
X^^' are operators of rank p in spin space and coordinate 
space, respectively. A contribution from second-order in 
perturbation theory to a state (p can then be written as 



AE, 



(11) 



where ipj defines an intermediate state with energy de- 
nominator AEj. The summation is done over all inter- 
mediate states V'j • As far as tpj varies in this summation, 
within a fixed configuration with respect to harmonic- 
oscillator (HO) shells, AEj remains constant due to the 
degeneracy of single-particle energies in a given HO shell. 
We mention that the usage of non-degenerate perturba- 
tion theory yields only small changes. 

Such a configuration for a given HO shell is denoted 
by S. As AEj is a constant within a fixed shell S, we 
label it as AEs ■ Note that S corresponds to a part of the 
Q-space, while (j) is in the P-space. The term 77(0) can 
then be decomposed into contributions from individual 
S"s as 



AEs 



(12) 



where 



C(<^,5) = E(,/)|iJi|V,)(V^,|Hi| 



(13) 



For a given S, all tpjs are included, and the summation 
can be replaced by the closure relation as 



C(0, S) ^ mH,H,}s\ 



(14) 



where the parentheses { }s are introduced to indicate 
that the second Hi changes (j> to an S'-configuration state 
in the Q-space and the first Hi moves it back to state cf> 
in the P-space. In other words, HiHi in this equation 
cannot be a simple product, but a certain contraction is 
needed as we shall show soon. 

By utilizing Eq. pH)) . we obtain 



{HiHijs = Wp.^pAiU'^'"^ ■ X^P^'^){U^P^-^ ■ X(^'^')}5 



Pl,P2 



fc=0.1.2 







P2 






Pi 


P2 








[ 










(0), 



k } { [U^P^^ X X [^^^'^ X X^P'')]^''^ }s , (15) 



where the terms in curly brackets are 9j symbols and 
k implies the rank of the recoupling. The operator 
X U'^P^^^'^^s acts in the P-space as a rank-A: two- 
body operator in spin space, while {[X^^^^ x X^^'^'J^'^^jg 



acts as a rank-A two-body operator in coordinate space. 
Because the contraction due to the elimination of the Q- 
spacc does not affect the angular momentum properties, 
the variable k = 0, 1,2 represents induced central, spin- 
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FIG. 4: (color online.) The tensor force monopole component of the effective interaction for the shell model obtained by the 
Q-box expansion to second and third order in the interaction, starting from Viowfc(A = 2.1 fm~^). The tensor-force component 
is obtained using the decomposition of Eq. ([7|. (a) T = forces in p/-shell, (b) T = 1 forces in p/-shell, (c) T = forces in 
sd-shell and (d) T = 1 forces in sd-shell. 



orbit and tensor forces in the P-space, respectively. 

Since we are mainly interested in the tensor compo- 
nent, we focus on the case of = 2, with the obvious 
restriction pi + P2 > 2. Since the above 9j symbols is 
proportional to 1/^ {2pi + l)(2p2 + 1), it is easy to con- 
vince oneself that the central force component receives 
the largest contribution from the 9j symbol. Further- 
more, for our analyses it is important to keep in mind 
that the expectation value of the central component is 
the largest in absolute value, the tensor component the 
second largest and the spin-orbit term gives rise to the 
smallest contribution to the renormalized V^owfe interac- 
tion. 

From these considerations, for fc = 2 the largest con- 
tribution comes from the combination pi = 0, p2 = 2 or 
Pi = 2, p2 = in Eq. ([TU|) . that is either a central-tensor 
or a tensor-central combination. Let us now discuss this 
case. We assume without loss of generality that the ten- 
sor component of Hi acts on the ket state of the matrix 
element being considered. While the central component 
of Hi acts afterward on this state, we can also consider 
that this central force acts to the left on the bra state. We 
then take the overlap between these two states by con- 
sidering one by the tensor on the ket side, and the other 
by the central on the bra side. These two states are sum- 
rule states for the two forces within the S'-configuration 



space. As the central force and the tensor force are very 
different in nature, such sum-rule states are very differ- 
ent from each other in general, leading to a very small 
overlap. This is the main reason why the combination 
of the central force and the tensor force produces small 
contributions. 

This argument does not hold for the case where the 
tensor component of Hi acts twice in the term to second- 
order in perturbation theory. However, due to the angu- 
lar momentum coupling, the product of two tensor forces 
(pi = p2 = 2 in Eq. (|10p ) yield small contributions to 
the k ~ 2 terms of Eq. ([TU|) . For higher orders, other 
tensor-force components may show up, but there is no 
mechanism to enhance their contributions. 

The small contribution of the tensor force in MBPT 
can be viewed to be reasonable also under the following 
intuitive picture: after multiple actions of the forces, the 
spin dependence is smeared out, and only the distance 
between two interacting nucleons becomes the primary 
factor to the whole processes. This results in the domi- 
nance of the induced effective interaction by the central 
components and yields only a minor change in the tensor 
component. 

It is instructive to study in more detail the contribu- 
tions to second-order in perturbation theory. To do so, 
we single out the by far largest second-order term, namely 
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the so-called core-polarization term, depicted as diagram 
(c) in Fig. [TJ For the core-polarization diagram we can 
show that the contribution to the tensor force vanishes by 



J 



simple angular momentum algebra arguments. The con- 
tribution to a specific core-polarization matrix element 
can then be written as 



= {amapmp\Vc\cmchmh){hmhbrni,\VT\p'mpdrnd) / A.E 

p.nipji.nih 

= ^ ^ ^ {amapmp\Vc\cmchm}i){hm}ibmb\VT\pmpdmd)^ I ^Es 

np,lp,nh,lh jp.mp,jh,mh 

= ^ ( ^ {amanplpmiprnsp\Vc\cmcnhlhrnihrnsh) 

np,lp,nh,lh rnip,m^p,vaih,m^h 



x{nhlhmhmshbrnb\VT\nplpmipmspdmd)j/AEs, 
I 



(16) 



where vp^~'^^ is the induced tensor force, Vr and Vc 
are the tensor force and central force components from 
Hi, respectively. With a harmonic oscillator basis, the 
term A.Es represents an energy denominator which is 
constant for a given set of quantum numbers np,nh-,lp 
and Ih, as discussed above. Here a = {na,la, ja)i and 
nia denotes magnetic substate of la- Note that the two- 
body states are not antisymmetrized. The states p and 
h represent particle and hole states, respectively. In the 
third line of the equation, only particle and hole states 
arc transformed to the Is coupling scheme. Note that 
the intermediate states are summed up to fulfill spin- 
saturation within each HO major shell. 

We can divide the contribution into two different types 
according to the spin dependence of the central force. 
One comes from the terms whose central force part Vc 
includes a ■ a (type I) and the other does not (type H). 
With our summation tailored to a spin-saturated core or 
an excluded Q-space with all spin-orbit partners, we can 
prove that type H contributions always vanish because 
the first factor is diagonal with respect to spin, that 
is rrisp = rrish, and the second factor is zero when we 
sum over spin-saturated contributions. Therefore, only 
a spin-dependent central force results in non-vanishing 
contributions to the tensor force for higher-order terms 
in V^f^. Finally, the contribution to the tensor force 
from the spin-dependent central force is quite small be- 
cause the spin-dependent central force is generally by far 
smaller than the spin-independent central force in mod- 
ern realistic NN potentials. 

In conclusion, medium effects produce minor contribu- 
tions to the tensor-force component, resulting in a tensor- 
force component that is dominated by the bare NN in- 
teraction. Our hypothesis about Renormalization Persis- 
tency is fulfilled to a good extent by the tensor force. 

Finally, although our analyses has been performed 
within one major shell only, one should note that this 
persistency of the tensor-force component via a MBPT 
renormalization should also hold in for model spaces 



which span several shells. 



VI. TWO-BODY MATRIX ELEMENTS OF THE 
TENSOR FORCE 

In this section we study further renormalization prop- 
erties of the tensor force by including higher multipole 
components. 

Figure [5] shows the diagonal and non-diagonal matrix 
elements of the bare tensor force from the AV8' poten- 
tial, the renormalized Viowfe interaction (A = 2.1 fm^^) 
and VJI^ obtained by the Q-box expansion up to the 
third order with folded diagrams to infinite order start- 
ing from AV8' interaction. This is similar to what was 
done in Figs. [5] and |31 In this figure, panels (a) to (d) 
stand for the sd-shell matrix elements. The diagonal ma- 
trix elements are shown in panels (a) and (b) while the 
non-diagonal elements are shown in panels (c) and (d). 
Note that the diagonal matrix elements {jajb\V\jajb) jt 
are specified by the quantum numbers ja, jt, and twice 
the total angular momentum J and total isospin T. The 
non-diagonal matrix elements {jajb\V\jcjd) jt are speci- 
fied by jai jb, jci id, J and T. The corresponding num- 
bers for the p/-shell are shown in panels (e) through (j). 
In both the sd-shell and the p/-shell, the patterns are the 
same for all approaches to the effective interactions and 
thus the R-Persistency is approximately fulfilled. In par- 
ticular, for the Viowfc renormalization procedure, we can 
hardly sec any difference between the bare tensor force 
and the tensor force in the effective interaction Mowfc ■ For 
the diagonal matrix elements, we can see small differ- 
ences between the final Q-box and the bare tensor force, 
however, it does not contradict the results with respect 
to the monopole component discussed above, mainly be- 
cause only matrix elements with small values of the total 
angular momentum display sizable differences. Since the 
monopole terms are weighted by 2 J -I- 1, matrix elements 
with larger values of the total angular momentum J carry 
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a much larger weight in Eq. ([TJ. In non-diagonal matrix 
elements, we see somewhat larger differences. Their role 
in shell-model calculations needs to be investigated fur- 
ther. A spin-tensor analysis along these lines was made 
recently by Smirnova et al [Toj . 



VII. ANALYSIS OF OTHER INTERACTION 
MODELS 

In the previous sections, we calculated effective inter- 
actions starting from the AV8' interaction, using a renor- 
malized interaction and many-body perturbation theory. 
We found that the R-Persistency of the tensor force holds 
for all these renormalization procedures. An obvious 
question is whether or not the R-Pcrsistency holds for 
other interaction models as well. In this section we ad- 
dress this issue as well. 

We employ here another frequently used realistic in- 
teraction, xN'^LO, as an example The xN'^LO 
interaction has a relatively smaller coupling between 
low-momentum and high-momentum modes compared 
with the AV8' potential. In Fig. [S] we show the 
monopole part of the tensor force of the xN'^LO bare 
potential and VlowfcS with several cutoff parameters 
A(1.0 fm"\ 2.1 fm"^ and 5.0 fm"^). These resuhs should 
be compared with the corresponding ones obtained with 
the AV8' interaction shown in Fig.[2l Figure [7] shows the 
tensor- force monopole of Vjj'^, corresponding to Fig. 0] 
for the AV8' interaction, starting from the xN^LO inter- 
action. Finally, Fig. |5] shows the multipole components 
of effective interactions, corresponding to Fig. [5] for the 
AV8' interaction. 

In all the figures, we can conclude that all the fea- 
tures we discussed for the AV8' interaction pertain to 
the xN^LO interaction model as well. 



bation theory in order to obtain an effective shell-model 
interaction. This has lead us to introduce the concept of 
renormalization persistency (R-Persistency) in the study 
of effective interactions. We studied the R-persistency 
of both renormalization procedures and showed via nu- 
merical studies, their intuitive general explanations and 
a detailed algebraic analysis of core-polarization terms 
in perturbation theory, that this is a very robust process. 
We have also shown that the R-Persistency holds for two- 
body matrix elements including higher multipole compo- 
nents of the tensor force, although the deviation increases 
somewhat if multipole components are included in the 
comparison. We conclude that the two renormalization 
steps (one for short-range correlation and the other for in- 
medium effects) do not affect much cither the monopole 
nor the multipole components of the tensor force, apart 
from slight differences between them. Results obtained 
with two different interactions (AV8' and xN'^LO) lead us 
to the same conclusion, suggesting that the R-Persistency 
of the tensor force for low-momentum states is a robust 
feature. This applies also to other interaction models 
than those studied here. 

The short-range part of the tensor force enters the 
renormalization of the central force, in particular in the 
T = channel, producing on average an increased attrac- 
tion. Since the modification of the tensor force appears to 
be small, the central force carries most of the renormal- 
ization effects beyond first order in perturbation theory. 

Because the R-Pcrsistency of the tensor force in effec- 
tive interactions is a robust feature, it may give a sim- 
ple and concrete starting point for examining and con- 
structing effective interactions, especially phenomenolog- 
ical ones In particular, since the tensor force plays a 
significant role in the shell evolution for nuclear systems 
with either large neutron/proton or proton/neutron ra- 
tios, one can extract simple physics messages from com- 
plicated many-body systems. 



VIII. CONCLUSION 



In this work we have presented a detailed analysis 
of various contributions to the nuclear tensor force as 
function of different renormalization procedures, start- 
ing with state-of-the-art nucleon-nucleon (NN) interac- 
tions and ending up with effective interactions for the nu- 
clear shell model. The monopole part of the tensor force 
is weakly or barely affected by various renormalization 
procedures, which in our case are represented by a renor- 
malization of the bare interaction and many-body pertur- 



Acknowledgments 

We are very grateful to Professors R. Okamoto and 
H. Feldmeier for valuable discussions. This work is sup- 
ported in part by Grant-in-Aid for Scientific Research 
(A) 20244022 and also by Grant-in-Aid for JSPS Fellows 
(No. 228635), and by the JSPS Core to Core program 
"International Research Network for Exotic Femto Sys- 
tems" (EFES). 



[1] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. 

Rev. C 51, 38 (1995). 
[2] R. B. Wiringa and S. C. Pieper, Phys. Rev. Lett. 89, 

182501 (2002). 

[3] T. Otsuka, T. Suzuki, R. Fujimoto, H. Grawe, and 
Y. Akaishi, Phys. Rev. Lett. 95, 232502 (2005). 



[4] T. Otsuka, T. Suzuki, M. Honma, Y. Utsuno, N. Tsun- 
oda, K. Tsukiyama, and M. Hjorth- Jensen, Phys. Rev. 
Lett. 104, 012501 (2010). 

[5] R. Bansal and J. French, Physics Letters 11, 145 (1964). 

[6] A. Poves and A. Zuker, Physics Reports 70, 235 (1981). 

[7] S. Bogner, T. T. S. Kuo, L. Coraggio, A. Covello, and 



10 



N. Itaco, Phys. Rev. C 65, 051301 (2002). 
[8] S. K. Bogner, T. T. S. Kuo, and A. Schwenk, Physics 

Reports 386, 1 (2003). 
[9] T. Kuo and E. Osnes, Lecture Notes in Physics 364, 1 

(1990). 

[10] M. Hjorth- Jensen, T. T. S. Kuo, and E. Osnes, Physics 

Reports 261, 125 (1995). 
[11] E. Epelbaum, H.-W. Hammer, and U.-G. Meifiner, Rev. 

Mod. Phys. 81, 1773 (2009). 
[12] M. W. Kirson, Physics Letters B 47, 110 (1973). 
[13] E. Osnes and D. Strottman, Phys. Rev. C 45, 662 (1992). 
[14] B. A. Brown, W. A. Richter, R. E. Julies, and B. H. 

Wildenthal, Annals of Physics 182, 191 (1988). 
[15] A. Nogga, S. K. Bogner, and A. Schwenk, Phys. Rev. C 



70, 061002 (2004). 
[16] S. K. Bogner, A. Schwenk, T. T. S. Kuo, and G. E. 

Brown, arXiv nucl-th (2001). 
[17] A. Umeya and K. Muto, Phys. Rev. C 74, 034330 (2006). 
[18] R. Tamagaki and W. Watari, Progress of Theoretical 

Physics Supplement 39, 23 (1967). 
[19] N. A. Smirnova, B. Bally, K. Heyde, F. Nowacki, and 

K. Sieja, Physics Letters B 686, 109 (2010). 
[20] For the j = j' case, the definition is slightly different. 
[21] A short version of the present results was included in [j, 

while more substantial, deeper and wider discussions are 

given in this paper. 



5' O 



CD , , 

<H O 

00 n 

g 

o B- 



ti 
o 
ti 

I 

& 

O 

ti 



cab|Vj|cd> (MeV) 



<ab|Vj|cd> (MeV) 



«ab|VT|ab> (MeV) 



cab|Vj|cd> (MeV) 



o 
B 
o 



o 
B 

to 

55 

CD 
O 
c+ 

^' 



cab|Vj|ab> (MeV) 



-L o o ro 
ai-^aioai-^biroai 




la lb Ic Id 

1777 8- 

1757 6- 

1757 8- 

1755 8- 

1737 6- 
1 737 

1735 6- 

1735 8- 

15 7 7 4- 

1557 4- 

1557 6- 

15 5 5 4- 

1537 4- 

15 3 7 6- 

15 3 5 4- 

1535 6- 

1533 4- 
15 17 

13 7 7 4- 

13 5 7 2- 

13 5 7 4- 

13 5 5 4- 

13 3 7 4- 

13 3 5 2- 

13 3 5 4- 

13 3 3 4- 

13 15 4- 

117 7 0- 

115 5 0- 
11330- 



o 
o 

— t 

o 

CD 

cn 



o 

CD 



pf-shell 



sd-shell 



> 

03 



> 
03 



0.5 



-0.5 



0.5 



(a) T=0 forces 




bare (N3L0] 
A= 1.0fm" 
A = 2.1 fm" 
A = 5.0fm"^ 





CO 

m Q. 


o o <^ 
V- V- Q. Q. Q. 

Q. Q. t t !C 


Q. 
1 


_ (b)T=1 


forces 


bare(N3L0] 
A= I.Ofm" - 


□ 






A = 2.1 fm" - 








A = 5.0fm"^ 








in 


CO 


CO 








"a 


"a 


"a 


03 


03 


03 


ih 


CO 


in 




in 


CO 


"a 


T3 


T3 


0) 







1 

0.5 




9- ^ !G "S. Q. 2. 
iii *Q- R t t t 



in 
in 



iG !G Q. 



(d) T=1 forces 


bare (AV8') □ 

A = I.Ofm"] 

A = 2.1 fm" 

A = 5.0fm"'' 









in 


CO 


CO 










"a 


"a 


03 


03 


03 


in 


CO 


in 




in 


CO 


T3 


T3 


T3 


0) 







FIG. 6: (color online.) Tensor-force monopole of Viowfe starting from the xN^LO interaction with the same notation as in Fig 
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